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ABSTRACT 


A  computer  program  was  written  that  enables  the  design  of  parametric 
sonars.  This  program  accepts  as  inputs  temperature,  salinity,  depth,  estimates 
of  projector  area,  desired  secondary  source  level  and  secondary  frequency. 

The  program  computes  various  parametric  sonar  quantities  among  them  primary 
source  level,  directivity  index  and  primary  operating  frequency.  The  program 
actually  generates  a  matrix  of  possible  design  values  that  permit  the  designer 
to  choose  those  which  best  suit  his  needs  based  on  other  system  considerations. 

The  design  program  is  written  in  Fortran  V  for  use  on  the  Univac  1108. 

The  program  is  completely  general  and  any  of  the  input  parameters  can  be  varied 
while  holding  the  others  constant.  A  discussion  on  how  to  use  the  program  as 
well  as  a  sample  example  is  included.^ _ 
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INTRODUCTION 

From  the  viewpoint  of  the  individual  who  is  faced  with  the  design  of  a 
parametric  sonar,  the  calculations  involved  seem  repetitive  and  in  some  cases 
endless.  The  Mellen  and  Moffett*  curves  together  with  the  appropriate  equations 
given  in  the  reference  contain  all  of  the  necessary  information.  However,  the 
information  is  presented  in  such  a  way  as  to  make  it  easy  to  work  through  the 
curves  and  calculations  to  analyze  the  parametric  operation  of  existing  projec¬ 
tors  and  systems  but  it  is  difficult  and  not  straightforward  to  work  backward 
through  the  curves  and  equations  to  design  a  projector  system. 

2 

There  is  a  solution  which  is  offered  by  Moffett  that  uses  a  "load  line" 
type  of  technique  similar  to  that  used  in  vacuum  tube  circuit  design.  This  is 
good  for  a  small  number  of  possible  designs  of  a  given  parametric  sonar.  The 
method  requires,  however,  repetitive  computations  to  arrive  at  the  dimensionless 
parameters  l/(2)(AL)(RO)  for  each  possible  parametric  stepdown  ratio  (the  ratio 
of  the  parametric  difference  frequency  to  the  mean  projector  driving  frequency 
FO/F).  The  term  (AL)  is  absorption  in  nepers  per  meter  while  RO  is  the  Rayleigh 
distance.  Appendix  B  contains  a  complete  glossary  of  terms.  The  "load  line" 
method  is  presently  limited  by  the  number  of  families  of  curves  available  for 
the  different  stepdown  ratios  and  the  accuracy  of  interpolating  between  the 
given  curves  of  a  given  family. 

A  means,  therefore,  was  devised  where  the  whole  design  process  was  automated 
using  the  Univac  1108  computer.  In  essence,  computer  aided  design.  The  solution 
allows  the  designer  to  work  from  a  known  secondary  source  level  (LSS)  and  a 
known  secondary  frequency  (F)  for  a  range  of  values  of  projector  size  (A), 
primary  source  level  per  tone  (LSP)  and  a  given  stepdown  ratio  (FO/F).  The 
computer  program  will  build  a  matrix  of  possible  designs  that  can  then  be  compared 
with  other  factors  to  achieve  a  workable  and  realistic  design. 


BACKGROUND 

In  parametric  sonar  calculations,  two  distinct  and  different  problems 
arise.  One  is  that  of  the  analysis  of  existing  sonar  to  predict  their  parametric 
operating  characteristics.  The  other  is  the  design  of  parametric  sonars  having 
a  given  set  or  range  of  output  source  levels  and  frequencies. 

In  the  first  problem,  one  usually  knows  the  primary  operating  frequency 
(F0)|  the  primary  source  level  (LSP),  and  the  projector  area  (A).  From  these 
one  obtains  the  secondary  source  level  (LSS)  and  secondary  directivity  index 
(NDIS)  for  a  given  downshift  ratio  (FO/F)  by  using  the  Mellen  and  Moffett 
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curves  and  the  appropriate  formulae.  In  summary,  we  know 

FO,  LSP  and  A. 

We  find 

LSS,  FO/F,  NDIS. 

The  Mellen  and  Moffett  curves  and  the  associate  equations  readily  lend  them¬ 
selves  to  solving  this  problem  because  of  the  way  the  equations  and  the  curves 
are  set  up. 

The  second  problem  is  one  of  designing  a  parametric  sonar  starting  from 
"scratch"  where  one  only  knows  the  desired,  or  the  range  of  desired,  LSS,  F, 
and  NDIS  and  wants  to  find  FO,  LSP  and  A.  At  first  glance  one  would  say, 

"Why  not  just  work  backwards  through  the  equations  with  the  aid  of  the  Mellen 
and  Moffett  curves?"  Alas,  life  is  not  so  simple.  The  equations  depend  on 
a  knowledge  of  FO,  and  A.  In  other  words,  something  must  be  known  about  the 
projector  before  starting.  Unfortunately,  determining  FO,  A,  and  LSP  is  the 
goal  of  the  design  process.  This  is  just  the  opposite  of  the  analysis  previously 
discussed.  There  is  a  method  that  has  been  proposed  by  Moffett  that  utilizes 
the  "load  line"  technique  previously  mentioned.  This  method  is  excellent  when  an 
exact  FO  and  A  are  sought  for  a  given  LSS,  F,  and  NDIS.  The  method  becomes  time 
consuming  and  requires  tedious  repetitive  calculations  when  a  range  of  values 
is  sought  and  when  one  needs  numerous  possibilities  in  order  to  examine  and 
choose  an  optimum  solution  based  on  factors  other  than  just  parametric  sonar 
considerations.  What  is  needed  is  a  method  of  constructing  a  matrix  of  possible 
parametric  sonar  designs  for  the  designer  to  weigh  in  consort  with  associated 
system  parameters.  In  summary  for  this  situation  we  know 

LSS,  FO/F,  NDIS 

and  we  want  to  find 

FO,  LSP  and  A. 


SOLUTION 

The  solution  to  the  problem  is  computer  aided  design.  A  program  was 
written  that  allows  the  designer  to  vary  F,  A,  LSS,  and  FO/F  in  order  to 
construct  the  desired  design  matrix.  The  program  compilation  is  given  in 
Appendix  A.  This  program  is  versatile  enough  so  that  three  other  parameters 
temperature  (T),  salinity  (S)  and  depth  (D)  can  be  varied  in  coarse  steps  and 
their  effects  on  the  design  studied.  The  results  are  tabulated  and  two  on-line 
plots  are  possible.  The  results  of  a  sample  example  are  shown  in  Appendix  B. 
The  on-line  plots  can  be  of  any  two  variables  and  each  plot  can  be  altered  by 
changing  a  computer  card. 
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At  present,  one  plot  is  acoustic  power  in  dB  (PADB)  vs.  FO/F  for  a  given 
LSS  with  the  parameters  T,  S,  D,  F  and  A  held  constant.  Then  either  LSS,  F, 
or  A  can  be  changed  and  another  plot  made.  Thus,  one  can  examine  the  range 
of  possible  designs  that  are  within  the  desired  power  budget  and  select  a  rea¬ 
sonable  one.  The  second  plot  is  secondary  directivity  index  (NDIS)  vs  FO/F  for 
the  same  given  conditions  as  in  the  previous  plot.  From  this  the  designer  can 
select  the  necessary  quantities  for  a  desired  range  of  NDIS.  Normally,  many 
plots  will  be  produced  resulting  in  families  of  LSS  curves  with  NDIS  and  PADB 
plotted  against  FO/F  with  T,  S,  D,  constant  for  many  combinations  of  F  and  A. 

Once  a  set  of  parameters  is  decided  upon,  the  appropriate  exact  constants  can 
be  obtained  from  the  tabulation. 

One  thing  that  has  been  done  to  aid  in  plot  comparisons  is  to  force  the 
plots  to  a  convenient  common  scale.  This  was  done  by  the  use  of  two  dummy  points 
on  each  plot.  This  was  necessary  because  the  routine  as  originally  compiled  by 
Gordon^  automatically  scaled  the  axis  for  the  plotting  range.  For  the  desired 
comparisons  of  plots,  such  scaling  is  undesirable. 

The  program  is  outlined  in  a  simplified  flow  chart  shown  in  Figure  1.  It 
operates  as  follows!  The  inner  loop  computes  parametric  sonar  design  constants 
for  each  of  a  sequence  of  FO/F  values.  This  is  done  for  each  of  a  sequence  of 
LSS  values  in  the  input  data  (LSS1).  Next,  the  two  inner  loops  are  repeated  for 
a  sequence  of  values  for  A  and  finally  these  three  innermost  loops  are  repeated 
for  each  value  of  F  in  the  input  data.  These  four  loops  generate  a  matrix  of 
possible  designs  for  the  ranges  of  FO/F,  LSS,  A,  and  F  chosen.  Each  matrix  is 
built  up  for  constant  values  of  T,  S,  and  D.  The  values  for  T,  S,  and  D  can  be 
altered  by  changing  them  when  the  data  is  programmed  into  the  computer. 

The  plots  as  presently  compiled  plot  after  each  sequence  of  FO/F  for  a 
given  LSS.  Thus,  a  family  of  curves  of  different  LSS  values  is  generated. 

These  are  for  each  combination  of  T,  S,  D,  F,  and  A  and  are  plotted  with  FO/F 
as  the  horizontal  axis  on  each  plot.  The  vertical  axis  on  plot  number  1  is 
PADB  while  the  vertical  axis  on  plot  number  2  is  NDIS. 

A  more  detailed  flow  chart  is  shown  in  Figure  2,  It  shows  an  expansion  of 
the  computational  block  diagram  of  Figure  1.  Thus,  the  location  of  the  various 
calculations  are  shown  along  with  the  appropriate  tests  required  to  keep  the 
program  bounded.  Once  the  data  is  entered,  the  calculations  leading  to  the 
quantity  l/(2)(AL)(RO)  are  made  where  the  attenuation  loss  is  (AL)  and  the 
Rayleigh  distance  is  (RO).  This  quantity  l/(2)(AL)(R0)  together  with  the  FO/F 
and  a  quantity  X  is  entered  into  a  numerical  integration  subroutine  devised  by 
Goldstein  .  The  X  is  a  parameter  that  ties  the  integration  to  a  scaled  source 
level  (L*)  which  is  a  normalized  parametric  quantity  in  the  Mellen  and  Moffett 
theory.  The  output  of  the  numerical  integration  enters  into  several  simple 
computations,  the  results  of  which  are  tested  to  see  if  they  fall  within  the 
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proper  programmed  bounds.  If  the  tests  are  failed,  a  new  value  of  X  is  chosen 
and  the  integration  routine  is  redone  and  retested.  Depending  on  how  the  tests 
are  passed,  the  program  either  proceeds  to  calculate  further  parametric  sonar 
quantities  for  the  given  solution  of  the  numerical  integration  or  the  program 
recognizes  that  the  numerical  integration  has  searched  as  far  as  it  can.  In 
any  event,  the  program  will  proceed  to  readout  the  results  in  a  table  then 
recycle  to  the  next  FO/F  in  the  innermost  loop.  Once  the  desired  LSS  values  has 
been  completely  investigated,  the  computer  constructs  the  two  on-line  plots 
previously  mentioned.  The  program  then  recycles  until  all  possible  values  of 
F,  A,  LSS  and  FO/F  have  been  investigated  and  all  plots  completed.  The  program 
then  terminates.  The  detailed  flow  chart  (Figure  2)  references  equations  which 
are  tabulated  in  Table  I. 

In  essence,  the  program  takes  some  known  values  for  a  given  condition  and 
hunts,  by  means  of  a  numerical  integration  routine  and  specific  tests,  for  other 
needed  values  to  completely  describe  a  parametric  sonar.  Since  usually  there 
is  a  range  of  desired  values,  the  program  builds  a  matrix  of  possible  solutions. 
The  accuracy  of  these  solutions  depends  on  the  accuracy  of  the  parameter  X  used 
in  the  numerical  integration.  Presently,  the  solution  calculates  an  LSS  which 
is  compared  with  the  input  LSS  (LSSl).  The  calculated  value  has  a  tolerance  of 
+0.82  dB.  The  resultant  LSP,  parametric  gain  (G),  acoustic  power  (PA  and  PADB), 
and  primary  frequency  directivity  index  (NDIP)  all  have  a  tolerance  of  +0.41  dB. 
The  NDIS  has  a  +0.82  dB  tolerance. 


PROGRAM  OPTIONS 

The  design  program  has  certain  options  as  a  result  of  the  general  form  in 
which  it  is  written.  The  program  contains  four  nested  loops  any  of  which  can 
be  varied  or  held  constant  by  appropriate  input  data  on  the  input  data  cards. 

The  plots  can  be  varied,  however,  this  may  involve  repositioning  the  plot  in 
the  program  as  well  as  changing  two  program  cards.  The  user  may  have  to  redimen¬ 
sion  the  storage  associated  with  the  loops  preceding  the  plot  in  order  to  be  sure 
the  data  computed  is  retained  until  the  plot  is  called. 
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Equation  No.  Equation 


i.  X(I+1)  =  1*1XI  F0R  86  VALUES  FR0M  x  =  0.090909 


2^. 

FT  =  21.9  x  io6"(1520/(T+273))  kHz 

35. 

AL  =  (1/8.68) (1/914.4)  j 

(1.86x10-2)(S)(FT)(F0)2/  [(FT)2+(F0)2]J  + 

^2.68xlO“’2(FO)2/FT  ]  + 

0.1(F0)2/(1+(F0)2)1  J (1-6. 33233x1 05D) 

nepers/meter 

5 

4  . 

C  =  1449.2  +  4.623T  -  0.0546(1^)  +  1.391  (S-35)  +  0.017D  METERS/SECOND 

56. 

NDIP  =  10  L0G10(4irA  (F0)2(103)2/C2)  DB 

6. 

PADB  =  LSP  -  70.8  -  NDIP 

DB 

7. 

PA  =  ANTILOG  £(  l/lO )( LSP- 

-70.8-NDIP)J  WATTS 

81. 

NDIS  =  (NDIP)  +3-10  L0G10  [l  +  ((F0)/F)(2 *  (AL)(R0)  +  X)J  DB 

TABLE  I 


PROGRAM  USAGE  AND  SAMPLE  EXAMPLE 

In  order  to  use  the  program,  the  user  must  stack  appropriately  formatted 
data  cards  in  a  fixed  order  at  the  end  of  the  program.  There  are  six  different 
types  of  cards.  These  cards  will  now  be  discussed  in  order  from  front  to  back 
of  the  stack. 

The  first  type  of  card  contains  only  one  card  and  comes  first  in  the  data. 

It  is  formatted  into  4  fields  of  one  integer  number  per  field.  Each  number 
must  be  right  justified  in  a  field  width  of  five  (Fortran  V  statement  (15)). 

The  first  field  uses  columns  1  through  5  and  contains  the  number  of  F  values. 

The  second  field  uses  columns  6  through  10  and  contains  the  number  of  A  values. 
The  third  field  uses  columns  11  through  15  and  contains  the  number  of  LSS  values. 
The  fourth  field  uses  columns  16  through  20  and  contains  the  number  of  FO/F 
values  plus  1.  This  arrangement  of  fields  is  summarized  in  Table  2. 
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The  second  type  of  card  contains  only  one  card  and  it  is  the  2nd  card  in 
the  data.  It  is  formatted  into  3  fields.  Each  field  contains  a  number  that  is 
written  in  a  floating  point  format  which  is  right  justified  in  a  field  width  of 
10  with  a  4  decimal  place  accuracy  (Fortran  V  statement  (F10.4)).  The  first 
field  uses  columns  1  through  10  and  contains  the  value  for  T.  The  second  field 
uses  columns  11  through  20  and  contains  the  value  for  5.  The  third  field  uses 
columns  21  through  30  and  contains  the  value  for  D.  These  fields  are  also 
summarized  in  Table  2. 

The  third  through  sixth  type  of  cards  may  contain  more  than  1  card  for  each 
type  but  only  one  value  for  each  card.  Thus,  one  must  use  as  many  cards  for 
each  type  as  there  are  values  associated  with  that  type  and  the  cards  for  each 
type  must  be  grouped  together.  Each  number  is  written  in  a  floating  point 
format  which  is  right  justified  in  a  field  width  of  10  with  a  5  decimal  place 
accuracy  (Fortran  V  statement  (F10.5)).  Each  third  type  of  card  gives  a  value 
for  F.  Each  fourth  type  of  card  gives  a  value  for  A.  Each  fifth  type  of  card 
gives  an  input  value  for  LSS  (LSSl).  Lastly,  each  sixth  type  of  card  gives  a 
value  for  FO/F.  Each  of  these  field  layouts  are  summarized  in  Table  2. 


Card  Type 

Card  Columns 

Fortran  IV  Format 

Aqreement 

1 

1-5 

15 

No.  of  F  values 

6-10 

15 

No.  of  A  values 

11-15 

15 

No.  of  LSS  values 

16-20 

15 

No.  of  F0/F  values 

2 

1-10 

F10.4 

T  in  deg  C 

11-20 

F10.4 

S  in  PPT 

21-30 

F10.4 

D  in  meters 

3 

1-10 

F10.5 

F  in  kHz 

4 

1-10 

F10.5 

A  in  sq.  meters 

5 

1-10 

F10.5 

LSS  in  dB 

6 

1-10 

F10.5 

FO/F 

TABLE  2 


DATA  CARD  FORMATS 
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The  use  of  this  program  requires  the  input  of  data  from  a  program  stored 
on  tape  in  the  NUSC  New  London  Laboratory  Univac  1108  files.  This  is  tape 
U183.  Different  parameter  plots  may  be  made  by  simply  changing  the  call  to 
plot  (Call  Plot  A)  statements.  There  are  two  such  plots  in  the  program.  The 
plot  routine  can  be  eliminated  by  removing  the  two  call  to  plot  cards  which 
are  located  adjacent  to  each  other  in  the  program.  The  rest  of  the  program 
should  run  and  the  table  of  results  printed. 

A  sample  example  will  now  be  discussed.  Suppose  we  want  to  design  a 
parametric  sonar  that  has  the  following  specifications: 

LSS  =  90  dB//l|Jbar-meter 
F  =  3  kHz 

Ndi  =  30  dB  to  35  dB 

and  the  power  budget  is  such  that  we  wish  to  minimize  its  consumption.  Assume 
that  the  system  will  work  in  the  ocean  (S  =  35  ppt)  and  that  the  system  must  be 
capable  of  operating  in  the  winter  (T  =  7 °C)  on  the  surface 
(D  =  0).  The  data  is  programmed  as  shown  in  Figure  3.  The  tables  of  results 
are  shown  in  Appendix  B  along  with  3  sets  of  plots.  Examination  of  the  results 
shows  several  design  possibilities  all  within  the  region  of  a  dip  in  the  PADB 
plots.  If  it  were  not  possible  to  examine  so  large  a  quantity  of  points,  the  dip 
quite  possibly  would  go  unnoticed  because  there  is  a  tendency  for  the  unwitting 
designer  to  assumethat  increased  stepdown  ratio  means  increased  power  consump¬ 
tion.  Apparently,  this  is  not  always  true.  When  the  desired  points  are  isolated 
on  the  plots,  the  designer  then  can  go  to  the  tables  and  from  them  he  can  deter¬ 
mine  the  design  that  gives  the  desired  source  level  within  the  ND1S  restrictions. 
The  desired  design  for  the  sample  example  is  the  one  underlined  in  the  appro¬ 
priate  table  of  Appendix  B  and  encircled  on  each  of  the  associated  PADB  and  NDIS 
vs  FO/F  plots.  The  selected  design  has  the  following  parameters: 

FO/F.=  10 

FO  =  30  kHz 

LSP  =  131 . 2  dB//lpbar-meter 

ICIP  =  36.1  dB 

NDIS  =  34.4  dB 

and  PA  =  267.1  watts/each  primary  frequency. 

Other  related  quantities  can  be  obtained  from  the  data  tables.  For  different 
applications  these  quantities  may  assume  importance  and  thus  are  readily  avail¬ 
able  if  design  tradeoffs  become  necessary. 
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A  computer  aided  parametric  sonar  design  program  has  been  written  for  the 
UNIVAC  1108.  This  program  allows  the  designer  to  take  a  given  secondary  source 
level  (LSS),  secondary  frequency  (F)  and  secondary  directivity  index  (NDIS)  and 
compute  a  range  of  possible  parametric  sonar  designs  that  will  satisfy  his  needs. 
Thus,  the  selection  of  sonar  parameters  is  no  longer  limited  by  the  difficulty 
of  examining  a  range  of  possible  parametric  designs.  The  sonar  designer  can 
now  construct  a  matrix  of  possible  designs  then  base  the  final  selection  on 
which  of  these  designs  best  fits  the  other  systems  parameters  being  considered. 

By  means  of  computer  aided  design,  literally  hundreds  of  possible  designs  for 
a  given  situation  can  be  investigated  in  a  short  time. 
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APPENDIX  A 
PROGRAM  COMPILATION 


The  program  compilation  shown  here  is  complete  with  subroutines  except 
for  the  plot  subroutine.  That  particular  one  was  on  tape  and  was  not  compiled 
as  was  the  material  that  was  put  in  by  means  of  a  deck  of  cards. 
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SAMPLE  EXAMPLE  READOUT 


The  sample  example  tables  and  plots  are  shown  here.  These  are  as  they 
appear  on  line  out  of  the  high  speed  printer  associated  with  the  computer. 
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ABSTRACT 


A  computer  program  was  written  that  enables  the  design  of  parametric 
sonars.  This  program  accepts  as  inputs  temperature,  salinity,  depth,  estimates 
of  projector  area,  desired  secondary  source  level  and  secondary  frequency. 

The  program  computes  various  parametric  sonar  quantities  among  them  primary 
source  level,  directivity  index  and  primary  operating  frequency.  The  program 
actually  generates  a  matrix  of  possible  design  values  that  permit  the  designer 
to  choose  those  which  best  suit  his  needs  based  on  other  system  considerations. 

The  design  program  is  written  in  Fortran  V  for  use  on  the  Univac  1108. 

The  program  is  completely  general  and  any  of  the  input  parameters  can  be  varied 
while  holding  the  others  constant.  A  discussion  on  how  to  use  the  program  as 
well  as  a  sample  example  is  included.^ _  _ _ 

ADMINISTRATIVE  INFORMATION 

This  memorandum  was  prepared  under  Project  No.  A-614-19j,  Principal  Investi¬ 
gator,  Dr.  A.  J.  Van  Woerkum,  Code  TC.  - - - 

The  authors  of  this  memorandum  are  located  at  the  New  London  Laboratory, 
Naval  Underwater  Systems  Center,  New  London,  Connecticut  06320. 
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INTRODUCTION 

From  the  viewpoint  of  the  individual  who  is  faced  with  the  design  of  a 
parametric  sonar,  the  calculations  involved  seem  repetitive  and  in  some  cases 
endless.  The  Mellen  and  Moffett*  curves  together  with  the  appropriate  equations 
given  in  the  reference  contain  all  of  the  necessary  information.  However,  the 
information  is  presented  in  such  a  way  as  to  make  it  easy  to  work  through  the 
curves  and  calculations  to  analyze  the  parametric  operation  of  existing  projec¬ 
tors  and  systems  but  it  is  difficult  and  not  straightforward  to  work  backward 
through  the  curves  and  equations  to  design  a  projector  system. 

o 

There  is  a  solution  which  is  offered  by  Moffett  that  uses  a  "load  line" 
type  of  technique  similar  to  that  used  in  vacuum  tube  circuit  design.  This  is 
good  for  a  small  number  of  possible  designs  of  a  given  parametric  sonar.  The 
method  requires,  however,  repetitive  computations  to  arrive  at  the  dimensionless 
parameters  l/(2)(AL)(RO)  for  each  possible  parametric  stepdown  ratio  (the  ratio 
of  the  parametric  difference  frequency  to  the  mean  projector  driving  frequency 
FO/F).  The  term  (AL)  is  absorption  in  nepers  per  meter  while  RO  is  the  Rayleigh 
distance.  Appendix  B  contains  a  complete  glossary  of  terms.  The  "load  line" 
method  is  presently  limited  by  the  number  of  families  of  curves  available  for 
the  different  stepdown  ratios  and  the  accuracy  of  interpolating  between  the 
given  curves  of  a  given  family. 

A  means,  therefore,  was  devised  where  the  whole  design  process  was  automated 
using  the  Univac  1108  computer.  In  essence,  computer  aided  design.  The  solution 
allows  the  designer  to  work  from  a  known  secondary  source  level  (LSS)  and  a 
known  secondary  frequency  (F)  for  a  range  of  values  of  projector  size  (A), 
primary  source  level  per  tone  (LSP)  and  a  given  stepdown  ratio  (FO/F).  The 
computer  program  will  build  a  matrix  of  possible  designs  that  can  then  be  compared 
with  other  factors  to  achieve  a  workable  and  realistic  design. 


BACKGROUND 

In  parametric  sonar  calculations,  two  distinct  and  different  problems 
arise.  One  is  that  of  the  analysis  of  existing  sonar  to  predict  their  parametric 
operating  characteristics.  The  other  is  the  design  of  parametric  sonars  having 
a  given  set  or  range  of  output  source  levels  and  frequencies. 

In  the  first  problem,  one  usually  knows  the  primary  operating  frequency 
(FO),  the  primary  source  level  (LSP),  and  the  projector  area  (A).  From  these 
one  obtains  the  secondary  source  level  (LSS)  and  secondary  directivity  index 
(NDIS)  for  a  given  downshift  ratio  (FO/F)  by  using  the  Mellen  and  Moffett 
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curves  and  the  appropriate  formulae.  In  summary,  we  know 

FO,  LSP  and  A. 


We  find 


LSS,  FO/F,  NDIS. 


The  Mellen  and  Moffett  curves  and  the  associate  equations  readily  lend  them¬ 
selves  to  solving  this  problem  because  of  the  way  the  equations  and  the  curves 
are  set  up. 


The  second  problem  is  one  of  designing  a  parametric  sonar  starting  from 
"scratch”  where  one  only  knows  the  desired,  or  the  range  of  desired,  LSS,  F, 
and  NDIS  and  wants  to  find  FO,  LSP  and  A.  At  first  glance  one  would  say, 

"Why  not  just  work  backwards  through  the  equations  with  the  aid  of  the  Mellen 
and  Moffett  curves?"  Alas,  life  is  not  so  simple.  The  equations  depend  on 
a  knowledge  of  FO,  and  A.  In  other  words,  something  must  be  known  about  the 
projector  before  starting.  Unfortunately,  determining  FO,  A,  and  LSP  is  the 
goal  of  the  design  process.  This  is  just  the  opposite  of  the  analysis  previously 
discussed.  There  is  a  method  that  has  been  proposed  by  Moffett  that  utilizes 
the  "load  line"  technique  previously  mentioned.  This  method  is  excellent  when  an 
exact  FO  and  A  are  sought  for  a  given  LSS,  F,  and  NDIS.  The  method  becomes  time 
consuming  and  requires  tedious  repetitive  calculations  when  a  range  of  values 
is  sought  and  when  one  needs  numerous  possibilities  in  order  to  examine  and 
choose  an  optimum  solution  based  on  factors  other  than  just  parametric  sonar 
considerations.  What  is  needed  is  a  method  of  constructing  a  matrix  of  possible 
parametric  sonar  designs  for  the  designer  to  weigh  in  consort  with  associated 
system  parameters.  In  summary  for  this  situation  we  know 


LSS,  FO/F,  NDIS 


and  we  want  to  find 

FO,  LSP  and  A. 


SOLUTION 

The  solution  to  the  problem  is  computer  aided  design.  A  program  was 
written  that  allows  the  designer  to  vary  F,  A,  LSS,  and  FO/F  in  order  to 
construct  the  desired  design  matrix.  The  program  compilation  is  given  in 
Appendix  A.  This  program  is  versatile  enough  so  that  three  other  parameters 
temperature  (T),  salinity  (S)  and  depth  (D)  can  be  varied  in  coarse  steps  and 
their  effects  on  the  design  studied.  The  results  are  tabulated  and  two  on-line 
plots  are  possible.  The  results  of  a  sample  example  are  shown  in  Appendix  B. 
The  on-line  plots  can  be  of  any  two  variables  and  each  plot  can  be  altered  by 
changing  a  computer  card. 
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At  present,  one  plot  is  acoustic  power  in  dB  (PADB)  vs.  FO/F  for  a  given 
LSS  with  the  parameters  T,  S,  D,  F  and  A  held  constant.  Then  either  LSS,  F, 
or  A  can  be  changed  and  another  plot  made.  Thus,  one  can  examine  the  range 
of  possible  designs  that  are  within  the  desired  power  budget  and  select  a  rea¬ 
sonable  one.  The  second  plot  is  secondary  directivity  index  (NDIS)  vs  FO/F  for 
the  same  given  conditions  as  in  the  previous  plot.  From  this  the  designer  can 
select  the  necessary  quantities  for  a  desired  range  of  NDIS.  Normally,  many 
plots  will  be  produced  resulting  in  families  of  LSS  curves  with  NDIS  and  PADb 
plotted  against  FO/F  with  T,  S,  D,  constant  for  many  combinations  of  F  and  A. 

Once  a  set  of  parameters  is  decided  upon,  the  appropriate  exact  constants  can 
be  obtained  from  the  tabulation. 

One  thing  that  has  been  done  to  aid  in  plot  comparisons  is  to  force  the 
plots  to  a  convenient  common  scale.  This  was  done  by  the  use  of  two  dummy  points 
on  each  plot.  This  was  necessary  because  the  routine  as  originally  compiled  by 
Gordon"*  automatically  scaled  the  axis  for  the  plotting  range.  For  the  desired 
comparisons  of  plots,  such  scaling  is  undesirable. 

The  program  is  outlined  in  a  simplified  flow  chart  shown  in  Figure  1.  It 
operates  as  follows:  The  inner  loop  computes  parametric  sonar  design  constants 
for  each  of  a  sequence  of  FO/F  values.  This  is  done  for  each  of  a  sequence  of 
LSS  values  in  the  input  data  (LSSl).  Next,  the  two  inner  loops  are  repeated  for 
a  sequence  of  values  for  A  and  finally  these  three  innermost  loops  are  repeated 
for  each  value  of  F  in  the  input  data.  These  four  loops  generate  a  matrix  of 
possible  designs  for  the  ranges  of  FO/F,  LSS,  A,  and  F  chosen.  Each  matrix  is 
built  up  for  constant  values  of  T,  S,  and  D.  The  values  for  T,  S,  and  D  can  be 
altered  by  changing  them  when  the  data  is  programmed  into  the  computer. 

The  plots  as  presently  compiled  plot  after  each  sequence  of  FO/F  for  a 
given  LSS.  Thus,  a  family  of  curves  of  different  LSS  values  is  generated. 

These  are  for  each  combination  of  T,  S,  D,  F,  and  A  and  are  plotted  with  FO/F 
as  the  horizontal  axis  on  each  plot.  The  vertical  axis  on  plot  number  1  is 
PADB  while  the  vertical  axis  on  plot  number  2  is  NDIS. 


A  more  detailed  flow  chart  is  shown  in  Figure  2.  It  shows  an  expansion  of 
the  computational  block  diagram  of  Figure  1.  Thus,  the  location  of  the  various 
calculations  are  shown  along  with  the  appropriate  tests  required  to  keep  the 
program  bounded.  Once  the  data  is  entered,  the  calculations  leading  to  the 
quantity  l/(2)(AL)(RO)  are  made  where  the  attenuation  loss  is  (AL)  and  the 
Rayleigh  distance  is  (RO).  This  quantity  l/(2)(AL)(R0)  together  with  the  FO/F 
and  a  quantity  X  is  entered  into  a  numerical  integration  subroutine  devised  by 
Goldstein4.  The  X  is  a  parameter  that  ties  the  integration  to  a  scaled  source 
level  (L*)  which  is  a  normalized  parametric  quantity  in  the  Mellen  and  Moffett 
theory.  The  output  of  the  numerical  integration  enters  into  several  simple 
computations,  the  results  of  which  are  tested  to  see  if  they  fall  within  the 
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proper  programmed  bounds.  If  the  tests  are  failed,  a  new  value  of  X  is  chosen 
and  the  integration  routine  is  redone  and  retested.  Depending  on  how  the  tests 
are  passed,  the  program  either  proceeds  to  calculate  further  parametric  sonar 
quantities  for  the  given  solution  of  the  numerical  integration  or  the  program 
recognizes  that  the  numerical  integration  has  searched  as  far  as  it  can.  In 
any  event,  the  program  will  proceed  to  readout  the  results  in  a  table  then 
recycle  to  the  next  FO/F  in  the  innermost  loop.  Once  the  desired  LSS  values  has 
been  completely  investigated,  the  computer  constructs  the  two  on-line  plots 
previously  mentioned.  The  program  then  recycles  until  all  possible  values  of 
F,  A,  LSS  and  FO/F  have  been  investigated  and  all  plots  completed.  The  program 
then  terminates.  The  detailed  flow  chart  (Figure  2)  references  equations  which 
are  tabulated  in  Table  I. 

In  essence,  the  program  takes  some  known  values  for  a  given  condition  and 
hunts,  by  means  of  a  numerical  integration  routine  and  specific  tests,  for  other 
needed  values  to  completely  describe  a  parametric  sonar.  Since  usually  there 
is  a  range  of  desired  values,  the  program  builds  a  matrix  of  possible  solutions. 
The  accuracy  of  these  solutions  depends  on  the  accuracy  of  the  parameter  X  used 
in  the  numerical  integration.  Presently,  the  solution  calculates  an  LSS  which 
is  compared  with  the  input  LSS  (LSS1).  The  calculated  value  has  a  tolerance  of 
+0.82  dB.  The  resultant  LSP,  parametric  gain  (G),  acoustic  power  (PA  and  PADB), 
and  primary  frequency  directivity  index  (NDIP)  all  have  a  tolerance  of  +0.41  dB. 
The  NDIS  has  a  +0.82  dB  tolerance. 


PROGRAM  OPTIONS 

The  design  program  has  certain  options  as  a  result  of  the  general  form  in 
which  it  is  written.  The  program  contains  four  nested  loops  any  of  which  can 
be  varied  or  held  constant  by  appropriate  input  data  on  the  input  data  cards. 

The  plots  can  be  varied,  however,  this  may  involve  repositioning  the  plot  in 
the  program  as  well  as  changing  two  program  cards.  The  user  may  have  to  redimen¬ 
sion  the  storage  associated  with  the  loops  preceding  the  plot  in  order  to  be  sure 
the  data  computed  is  retained  until  the  plot  is  called. 
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Equation  No.  Equation 


1.  X(i+i)  =  l.lXj  FOR  86  VALUES  FROM  X  =  0.090909 


25. 

FT  =  21.9  x  i06-(1520/(T+273))  kHz 

35. 

AL  =  (1/8.68) (1/914.4)  j 
^2.68x10”2(FO)2/FT  J  + 
nepers/meter 

( 1 . 86x1 0"2 )  ( S )  ( FT)  ( FO) 2/  [( FT ) 2+(  FO) 2]  ]  + 

1  1  J 

0.1(F0)2/(1+(F0)2)1  [  (1-6. 33233x1 05D) 

5 

4  . 

C  =  1449.2  +  4.623T  -  0. 

0546(T2)  +  1.391  ( S-35 )  +  0.017D  METERS/ SECOND 

56. 

NDIP  =  10  L0G10(4tTA  (F0)2(103)2/C2)  DB 

6. 

PADB  =  LSP  -  70.8  -  NDIP 

DB 

7. 

PA  =  ANTILOG  [(1/10) (LSP 

-70.8-NDIP)]  WATTS 

81. 

NDIS  =  (NDIP)  +  3-10  L0G10  [l  +  ((F0)/F)(2ir  (AL)(R0)  +  X)J  DB 

TABLE  I 


PROGRAM  USAGE  AND  SAMPLE  EXAMPLE 

In  order  to  use  the  program,  the  user  must  stack  appropriately  formatted 
data  cards  in  a  fixed  order  at  the  end  of  the  program.  There  are  six  different 
types  of  cards.  These  cards  will  now  be  discussed  in  order  from  front  to  back 
of  the  stack. 

The  first  type  of  card  contains  only  one  card  and  comes  first  in  the  data. 

It  is  formatted  into  4  fields  of  one  integer  number  per  field.  Each  number 
must  be  right  justified  in  a  field  width  of  five  (Fortran  V  statement  (15)). 

The  first  field  uses  columns  1  through  5  and  contains  the  number  of  F  values. 

The  second  field  uses  columns  6  through  10  and  contains  the  number  of  A  values. 
The  third  field  uses  columns  11  through  15  and  contains  the  number  of  LSS  values. 
The  fourth  field  uses  columns  16  through  20  and  contains  the  number  of  FO/F 
values  plus  1.  This  arrangement  of  fields  is  summarized  in  Table  2. 


9 


TM  No. 
TD1X-33-73 


The  second  type  of  card  contains  only  one  card  and  it  is  the  2nd  card  in 
the  data.  It  is  formatted  into  3  fields.  Each  field  contains  a  number  that  is 
written  in  a  floating  point  format  which  is  right  justified  in  a  field  width  of 
10  with  a  4  decimal  place  accuracy  (Fortran  V  statement  (F10.4)).  The  first 
field  uses  columns  1  through  10  and  contains  the  value  for  T.  The  second  field 
uses  columns  11  through  20  and  contains  the  value  for  5.  The  third  field  uses 
columns  21  through  30  and  contains  the  value  for  D.  These  fields  are  also 
summarized  in  Table  2. 

The  third  through  sixth  type  of  cards  may  contain  more  than  1  card  for  each 
type  but  only  one  value  for  each  card.  Thus,  one  must  use  as  many  cards  for 
each  type  as  there  are  values  associated  with  that  type  and  the  cards  for  each 
type  must  be  grouped  together.  Each  number  is  written  in  a  floating  point 
format  which  is  right  justified  in  a  field  width  of  10  with  a  5  decimal  place 
accuracy  (Fortran  V  statement  (F10.5)).  Each  third  type  of  card  gives  a  value 
for  F.  Each  fourth  type  of  card  gives  a  value  for  A.  Each  fifth  type  of  card 
gives  an  input  value  for  LSS  (LSSl).  Lastly,  each  sixth  type  of  card  gives  a 
value  for  FO/F.  Each  of  these  field  layouts  are  summarized  in  Table  2. 


Card  TvDe 

Card  Columns 

Fortran  IV  Format 

Aqreement 

1 

1-5 

15 

No.  of  F  values 

6-10 

15 

No.  of  A  values 

11-15 

15 

No.  of  LSS  values 

16-20 

15 

No.  of  FO/F  values 

2 

1-10 

F10.4 

T  in  deg  C 

11-20 

F10.4 

S  in  PPT 

21-30 

F10.4 

D  in  meters 

3 

1-10 

F10.5 

F  in  kHz 

4 

1-10 

F10.5 

A  in  sq.  meters 

5 

1-10 

F10.5 

LSS  in  dB 

6 

1-10 

F10.5 

FO/F 

TABLE  2 


DATA  CARD  FORMATS 
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The  use  of  this  program  requires  the  input  of  data  from  a  program  stored 
on  tape  in  the  NUSC  New  London  Laboratory  Univac  1108  files.  This  is  tape 
U183.  Different  parameter  plots  may  be  made  by  simply  changing  the  call  to 
plot  (Call  Plot  A)  statements.  There  are  two  such  plots  in  the  program.  The 
plot  routine  can  be  eliminated  by  removing  the  two  call  to  plot  cards  which 
are  located  adjacent  to  each  other  in  the  program.  The  rest  of  the  program 
should  run  and  the  table  of  results  printed. 

A  sample  example  will  now  be  discussed.  Suppose  we  want  to  design  a 
parametric  sonar  that  has  the  following  specifications: 

LSS  =  90  dB//lpbar-meter 
F  =  3  kHz 

Ndi  =  30  dB  to  35  dB 

and  the  power  budget  is  such  that  we  wish  to  minimize  its  consumption.  Assume 
that  the  system  will  work  in  the  ocean  (S  =  35  ppt)  and  that  the  system  must  be 
capable  of  operating  in  the  winter  (T  =  7°C)  on  the  surface 
(D  =  0).  The  data  is  programmed  as  shown  in  Figure  3.  The  tables  of  results 
are  shown  in  Appendix  B  along  with  3  sets  of  plots.  Examination  of  the  results 
shows  several  design  possibilities  all  within  the  region  of  a  dip  in  the  PADB 
plots.  If  it  were  not  possible  to  examine  so  large  a  quantity  of  points,  the  dip 
quite  possibly  would  go  unnoticed  because  there  is  a  tendency  for  the  unwitting 
designer  to  assumethat  increased  stepdown  ratio  means  increased  power  consump¬ 
tion.  Apparently,  this  is  not  always  true.  When  the  desired  points  are  isolated 
on  the  plots,  the  designer  then  can  go  to  the  tables  and  from  them  he  can  deter¬ 
mine  the  design  that  gives  the  desired  source  level  within  the  NDIS  restrictions. 
The  desired  design  for  the  sample  example  is  the  one  underlined  in  the  appro¬ 
priate  table  of  Appendix  B  and  encircled  on  each  of  the  associated  PADB  and  NDIS 
vs  FO/F  plots.  The  selected  design  has  the  following  parameters: 

F0/F.=  10 

FO  =  30  kHz 

LSP  =  131 . 2  dB//lpbar-meter 

IPIP  =  36.1  dB 

NDIS  =  34.  4  dB 

and  PA  =  267.1  watts/each  primary  frequency. 

Other  related  quantities  can  be  obtained  from  the  data  tables.  For  different 
applications  these  quantities  may  assume  importance  and  thus  are  readily  avail¬ 
able  if  design  tradeoffs  become  necessary. 
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CONCLUSION 

A  computer  aided  parametric  sonar  design  program  has  been  written  for  the 
UNIVAC  1108.  This  program  allows  the  designer  to  take  a  given  secondary  source 
level  (LSS),  secondary  frequency  (F)  and  secondary  directivity  index  (NDIS)  and 
compute  a  range  of  possible  parametric  sonar  designs  that  will  satisfy  his  needs. 
Thus,  the  selection  of  sonar  parameters  is  no  longer  limited  by  the  difficulty 
of  examining  a  range  of  possible  parametric  designs.  The  sonar  designer  can 
now  construct  a  matrix  of  possible  designs  then  base  the  final  selection  on 
which  of  these  designs  best  fits  the  other  systems  parameters  being  considered. 

By  means  of  computer  aided  design,  literally  hundreds  of  possible  designs  for 
a  given  situation  can  be  investigated  in  a  short  time. 
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APPENDIX  A 
PROGRAM  COMPILATION 


The  program  compilation  shewn  here  is  complete  with  subroutines  except 
for  the  plot  subroutine.  That  particular  one  was  on  tape  and  was  not  compiled 
as  was  the  material  that  was  put  in  by  means  of  a  deck  of  cards. 
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appendix  b 

SAMPLE  EXAMPLE  READOUT 


The  sample  example  tables  and  plots  are  shown  here.  These  are  as  they 
appear  on  line  out  of  the  high  speed  printer  associated  with  the  computer. 
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